Novel phase-separation transition in one-dimensional driven models 
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A class of models of driven diffusive systems which is shown to exhibit phase separation in d = 1 
dimensions is introduced. Unlike all previously studied models exhibiting similar phenomena, here 
the phase separated state is fluctuating in the bulk of the macroscopic domains. The nature of 
the phase transition from the homogeneous to the phase separated state is discussed in view of a 
recently introduced criterion for phase separation in one-dimensional driven systems. 
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One dimensional driven diffusive systems have at- 
tracted considerable attention in recent years PJ. It 
has been demonstrated in numerous studies that unlike 
systems in thermal equilibrium, certain driven diffusive 
models with local noisy dynamics do exhibit phenomena 
like phase transitions and phase separation. More re- 
cently a general criterion for the existence of phase sep- 
aration in driven one dimensional models has been in- 
troduced (|] . The criterion relates the existence of phase 
separation in a given model to the rate at which domains 
of various sizes exchange particles with each other. As- 
suming that for a domain of length n this rate is given by 
the steady state current J n which flows through it, phase 
separation was suggested to exist only in the following 
cases: either the current vanishes in the thermodynamic 
limit. 



J n — > as n 



(Case A) 



(1) 



or the behavior of the current for large domains is of the 
form 



Jn ~ Joo{l + b/n a ) 



(Case B) 



(2) 



for either a < 1 and b > 0, or for a = 1 and b > 2. 

The nature of the phase separated states is rather dif- 
ferent in the two cases. In Case A the phase separated 
states were found to be of a rather simple nature, char- 
acterized by coexistence of pure domains, each consisting 
of a single type of particles. Thus, the particle density in 
the interior of a domain is non-fluctuating. Density fluc- 
tuations are limited to finite regions around the domain 
boundaries. Such steady states were termed strongly 
phase separated. Moreover, in this case phase separa- 
tion is expected to take place at any density, no matter 
how small. On the other hand, in Case B the phase sepa- 
rated state is expected to be fluctuating in the bulk of the 
macroscopic domains, as is normally expected in a noisy 
system. It exists only at high enough densities, while at 
low densities the system is homogeneous. This phase was 
termed condensed as the mechanism of the transition is 
similar to that of the Bose-Einstein condensation. 

So far, all one-dimensional models found to phase sep- 
arate are of type A |j, 0, || , and thus they exhibit strong 



phase separation at any density. In these models more 
than one species of particles is involved. In a recent study 
by Arndt et. al. (AHR) j| an interesting two species 
driven model was introduced. It was suggested, based 
on numerical simulations, that the model exhibits a con- 
densed phase separated state, whereby the particle den- 
sities fluctuate in the interior of the coexisting domains, 
and not just at the domain boundaries. In this state, 
a region with a high density of particles of both species 
coexists with a low density region. Moreover, the model 
has non-vanishing currents even in the thermodynamic 
limit. As in equilibrium phase separation it has been 
suggested that this state exists only at sufficiently high 
densities. However, a subsequent exact solution of the 
model H shows that what numerically seems like a con- 
densed state is in fact homogeneous, with a very large but 
finite correlation length. Further analysis of this model, 
in the light of the criterion suggested in [|j shows that the 
currents J n corresponding to this model are given by the 
form B, with er = 1 and b = 3/2 Therefore, according 
to the criterion, no phase separation takes place. 

Another example of a model which was suggested to 
exhibit phase separation into a fluctuating macroscopi- 
cally inhomogeneous state is the two-lane model intro- 
duced by Korniss et. al. R|. While numerical studies of 
the model indicate that such a phase exists in the model, 
studies of the current J n of finite domains suggests that it 
is of type B with a = 1 and b ~ 0.8 0, indicating, again, 
that no phase separation exists in this model. Thus the 
question of whether a phase separation of type B exists 
remains an intriguing open question. 

In this Letter we introduce a class of models which are 
demonstrated to be of type B, with a = 1 and b > 2. Ac- 
cording to the criterion conjectured in |2j this class is ex- 
pected to exhibit a phase transition to a phase separated 
condensed state. Thus at high densities these models 
exhibit a novel phase separation with non- vanishing cur- 
rents in the thermodynamic limit, and bulk fluctuations 
which are not restricted to the vicinity of the domain 
boundaries. To our knowledge, this is the first example 
of a genuine transition of this type in one-dimensional 
driven systems. 



We now define this class of models in detail. We con- 
sider a one-dimensional ring with L sites. Each site i 
can be either vacant (0) or occupied by a positive (+) 
or a negative (— ) particle (or charge). Positive parti- 
cles are driven to the right while negative particles are 
driven to the left. In addition to the hard-core repulsion, 
particles are subject to short-range interactions. These 
interactions are "ferromagnetic" , in the sense that parti- 
cles of the same kind attract each other. The dynamics 
conserves the number of particles of each species, N + 
and i\T_. The total density of particles in the system is 
p = (N + + N-)/L. The model is defined by a random- 
sequential local dynamics, whereby a pair of nearest- 
neighbor sites is selected at random, and the particles 
are exchanged with the following rates: 

H — -> — h with rate 1 - AH 

+ 0^0+ with rate a (3) 

— — > — with rate a . 

Here AH is the difference in the ferromagnetic interac- 
tions between the final and the initial configurations. We 
begin by considering a model with only nearest neighbor 
interactions, 

H = -e/4j2^s l+1 . (4) 

i 

Here Sj = +1 (—1) if site i is occupied by a + (— ) particle, 
and Si = if site i is vacant. The interaction parameter e 
satisfies < e < 1 to ensure positive transition rates. The 
model is a generalization of the Katz-Lebowitz-Spohn 
(KLS) model, introduced in |8| and studied in detail in 
p0[ , in which the lattice is fully occupied by charges and 
no vacancies exist. In this Letter we consider the case 
where the number of positive and negative particles is 
equal, N + = N— 

We will demonstrate that for a certain range of the 
parameters defining the dynamics, namely for e > 0.8 
and sufficiently large a (to be discussed below), a phase 
separation transition occurs as the density p is increased 
above a critical density p c . In the phase separated state 
a macroscopic domain, composed of positive and nega- 
tive particles, coexists with a fluid phase, which consists 
of small domains of particles (of both charges) separated 
by vacancies. Typical configurations obtained during the 
time evolution of the model starting from a random ini- 
tial configuration are given in Fig. @. This fi gure suggests 
that a coarsening process takes place, leading to a phase 
separated state as described above. However, this by it- 
self cannot be interpreted as a demonstration of phase 
separation in these models. The reason is that this be- 
havior may very well be a result of a very large but finite 
correlation length, as is the case in the AHR |) and the 
two- lane Q models discussed above ||. We thus apply 
the criterion pi in order to analyze the possible existence 
of phase separation in this model. 
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FIG. 1: Evolution of a random initial configuration of model 
(g) with nearest-neighbor interactions, on a ring of 200 sites. 
Here e = 0.9, a = 2, and the particle density is p — 0.5. 
Positive particles are colored black, and negative particles are 
colored grey. One hundred snapshots of the system are shown 
every 100 Monte-Carlo sweeps. 

To this end we note that a domain may be defined as an 
uninterrupted sequence of positive and negative particles 
bounded by vacancies from both ends. The current J n 
corresponding to such a domain of length n may thus be 
determined by studying an open chain, fully occupied by 
positive and negative particles, with entrance and exit 
rates a. This is just the one-dimensional KLS model on 
an open chain. Phase separation is expected to take place 
only for sufficiently large a. We consider a such that 
the system is in its maximal current state, whereby Joo 
assumes its maximum possible value, and is independent 
of a. 

To evaluate J n we first consider the KLS model on 
a ring of n sites with no vacancies. We then extend 
these results to study the behavior of an open chain. 
Since we are interested in the maximal current phase we 
consider equal number of positive and negative particles 
n + = n_ = n/2. It can be shown, as was done for the 
noisy Burger's equation jl^, |f3), that under quite general 
conditions, to be discussed below, the current J„ takes 
the following form for large n, 

Jn = Joo(l- ■ (5) 

\ 2,700 n ) 

Here A = d 2 J oa /dp 2 + is the second derivative of the cur- 
rent with respect to the density of positive particles p+ 
in the system. The compressibility analog k is defined by 
k = lhrin^oon -1 {{n 2 + ) — (n + ) 2 ), as calculated within a 
grand canonical ensemble, as explained below. This can 
be demonstrated by considering the current J n (n + ) for 
charge densities close to n + = n_ = n/2. Expanding 
J n (n+) in powers of An + = n + — n/2 one has 

Jn(n+) = J n (n/2) + J' n An+ + lj» (An + ) 2 (6) 

where the derivatives J' n and J„ are taken with respect 
to n + and evaluated at n/2. We average (|6|) over n + with 
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the steady state weights of a grand canonical ensemble. 
This is done by introducing a chemical potential p which 
ensures that the average density satisfies (n + ) = n/2. 
We find 



1 



(J„(n + ))„ = J n (n/2) + -J4' <(An+) 2 ) 



(7) 



Noting that (J n (n+)) /i is Joo in the n — > oo limit, and 
J n {n/2) is just J„, Eq. |^ is obtained. Here we made 
use of the fact that finite size corrections to (J„(n+)}^, 
resulting from the next to leading eigenvalue of the 
transfer-matrix of the steady-state distribution, are ex- 
ponentially small in n and may thus be neglected. The 
result of Eq. || is rather general, and is independent of 
the exact form of the steady-state particle distribution. 
This is provided that the weights of the microscopic con- 
figurations are local and thus the density and chemical 
potential ensembles are equivalent. 

In fact, an alternative way to derive (Q) is to consider 
the correspondence between the driven lattice-gas mod- 
els and the noisy Burger's equation or the Kardar-Parisi- 
Zhang (KPZ) equation for interface growth in 1 + 1 di- 
mensions plj . In these models J„ corresponds to the 
growth velocity of the interface. Eq. || has been derived 
in (l2[ where A is the coefficient of the non- linear 
term in the KPZ equation. The equivalence of the two 
alternative approaches relies on the fact that both k and 
A are invariant under renormalization transformations. 

The result (|J) can be used to evaluate J n for the KLS 
model. It has been shown || [n| that for a ring geometry 
the steady state weight of a configuration {t^} is 



P({n}) 



= e -m 



(8) 



with Tj = ±1 for positive and negative charges respec- 
tively, e 4/3 = (1 - e)/(l + e), and p serves as a chemical 
potential which controls the density of, say, the positive 
particles. The chemical potential p vanishes for the case 
n + = rt_. Using (||) expressions for «(e) and Joo(e) of 
this model have been obtained in |jTo|| . 

We now consider the KLS model in an open chain, 
which is the relevant geometry in applyingthe phase- 
separation criterion. It has been argued [jl2j that the 
finite size correction to the current of an open chain is 
given by the corresponding correction in a ring geometry, 
up to a universal multiplicative constant c which depends 
only on the boundary conditions. In the maximal current 
phase, c was found to be 3/2. Thus the current of an open 
system is given by (0) with a — 1 and 



6(e) 

Using the values of J OT 
we find 



A(e)/c(e) 
C 2J oc (e) 

and k obtained in [fL( 



(9) 

and c = 3/2 



6(e) 



3 (2 + e)v + 2e 
2 2(v + e) 



1 



1 



(10) 




FIG. 2: The coefficient 6(e), Eq. |l0| 

In figure || the coefficient 6(e) is plotted for < e < 1. 
This curve has been verified by direct numerical simula- 
tions of the KLS model on an open chain in the maximal 
current phase, demonstrating that the prefactor c indeed 
does not depend on e . Using ( |Toj ) it is readily seen that 
for e > 0.8 the value of 6 is larger than 2. 

According to the criterion conjectured in || one ex- 
pects phase separation to take place at high densities in 
model (||) for e > 0.8, as long as a is such that the KLS 
model is in the maximal current phase. This condensed 
phase separated state belongs to Case B of the criterion. 
We have carried out extensive numerical simulations of 
the dynamics of the model for various values of e. We 
find that for e < 0.4 no phase separation is observed. 
However, for e > 0.4 simulation of systems of sizes up 
to L = 10 6 show that the system evolves towards what 
seems to be a phase separated state at sufficiently large 
densities. We argue that a genuine phase separation takes 
place only for e > 0.8. On the other hand, the seemingly 
phase separation found in simulations for 0.4 < e < 0.8, 
is only a result of large but finite correlation lengths, as 
was found in the AHR and in the two-lane models. As 
pointed out in || such a behavior is related to correc- 
tions of order 1/n 2 and higher in the current (||). These 
corrections were shown to lead to a crossover with a very 
sharp increase in the correlation length, which could be 
erroneously interpreted as a genuine phase transition in 
numerical studies of finite systems. 

We now discuss the phase transition leading to the 
phase separated state. According to |2j the domain size 
distribution just below the transition takes the form 



V{n) 



(11) 



where £ is the correlation length, which diverges at the 
transition. The particle density in the system is related to 
£ by p/(l — p) =Y1 nr P( n )/J2 T ne critical density 

p c is given by this expression with £ — > oo. Note that 
with this form of the distribution function, p c is 1 in 
the limit 6 \ 2, and is a decreasing function of b. It is 



4 




FIG. 3: The coefficient b(8), as calculated from Monte-Carlo 
simulations of domains of sizes up to 1024. Data are shown 
for e = and e — 0.5 . 



straightforward to show that the divergence of the 
correlation length at the critical density is given by 



\P-Pcl~ , 2<6<3 



\P-Pc\~ 1 



(12) 



b > 3 



It is worthwhile noting that while /dp is continuous 
at the transition for 2 < b < 3, it exhibits a disconti- 
nuity for b > 3. The transition may thus be considered 
continuous for 2 < b < 3 and first-order for b > 3. 

In the model defined above b is found to satisfy 3/2 < 
b < 9/4. It is natural to ask whether larger values of b 
could be reached by increasing the range of the interac- 
tions. To answer this question we have extended model 
(Q) to include next-nearest-neighbor interactions as well, 
and consider 



H = -e/4]T 



V4£ 



SiS i+2 



(13) 



We have calculated the value of b as a function of S by 
Monte-Carlo simulations. This is done by measuring the 
current J n in an open system of size n, which is fully oc- 
cupied by positive and negative particles. At the bound- 
aries, the coupling to the rest of the system is modeled 
by injection of positive (negative) particles with rates a 
at the left (right). Simulating systems of size up to 1024 
enables us to fit the measured values of J n to the form 
(||) with cr=l, and to extract b. In figure |3| we plot b 
as a function of S, for e = and for e = 0.5. We find 
that by extending the range of the interactions one can 
increase b to values even larger than 3, where the phase 
separation transition is expected to be first order. 

In summary, a class of driven diffusive models in one- 
dimension is introduced and analyzed using a recently 
conjectured criterion for phase separation H. These 



models are shown to exhibit a novel type of phase sep- 
aration. In the phase separated state of these models 
the density is fluctuating in the bulk of the domains. 
Moreover, the models exhibit a homogeneous state at 
low densities, and a phase transition into the phase sep- 
arated state occurs at a critical density. The nature of 
the phase transition in these models is also discussed. 

While the validity of the criterion was proved for the 
AHR model, its general validity was conjectured based 
on some plausible assumptions on the behavior of the 
coarsening domains [^| . It would be of interest to analyze 
the class of models introduced in the present study by 
other analytical means, in order to verify the validity of 
the criterion. 
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